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ABSTRACT 

We report on the photospheric evolution of an intermediate-scale (« 4 Mm footpoint separation) 
magnetic bipole, from emergence to decay, observed in the quiet Sun at high spatial (O'.'S) and temporal 
(33 s) resolution. The observations were acquired by the IMaX imaging magnetograph during the first 
science flight of the Sunrise balloon-borne solar observatory. The bipole flux content is 6 x 10^^ Mx, 
representing a structure bridging the gap between granular scale bipoles and the smaller ephemeral 
regions. Footpoints separate at a speed of 3.5 km s~^ and reach a maximum distance of 4.5 Mm before 
the field dissolves. The evolution of the bipole is revealed to be very dynamic: we found a proper 
motion of the bipole axis and detected a change of the azimuth angle of 90° in 300 seconds. The overall 
morphology and behaviour are in agreement with previous analyses of bipolar structures emerging at 
granular scale, but we also found several similarities with emerging flux structures at larger scale. The 
flux growth rate is 2.6 x lO^'^ Mx s~^, while the mean decay rate is one order of magnitude smaller. We 
describe in some detail the decay phase of the bipole footpoints which includes break up into smaller 
structures, interaction with pre-existing flelds leading to cancellation but appears to be dominated by 
an as-yet unidentified diffusive process that removes most of the flux with an exponential flux decay 
curve. The diffusion constant (8 x lO'^ km^ s^^) associated with this decay is similar to the values 
used to describe the large scale diffusion in flux transport models. 
Subject headings: Sun: magnetic topology — Sun: photosphere — Sun: activity 



1. INTRODUCTION 

Magnetic fiux emergence in the solar photosphere 
involves a variety of spatial scales. It spans 
from the sr nall bipoles with fluxes of the order 
of 10^^ Mx (iMartinez Gonzalez fc Bellot RubTol 120091: 
iDanilovic et al.l |2010[) . that populate the quiet Sun at 
any stage of the solar activity cycle, to the large, com- 
plex active regions present dur ing solar maximu m, with 
fluxes of the order of 10^^ Mx (|Lites et al.l |1998). 

Bipolar structures with absolute flux content ranging 
from 10^* Mx to 5 x 10^^ Mx are usually referred to as 
ephemeral regions (ERs), due to the short lifetime of less 
than 24 ho urs that was estimated when they were first 
discovered (jHarvev fc Martinl ll973[ ). The opposite polar- 
ity components of ERs separate from each other during 
the growth period, reaching a distance of a few mega- 
meters in about half an hour with decreasing separation 
velocity ()Martin 1119881) . Then, ERs decay in a complex 
way, strongly dependent upon the surrounding magnetic 
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network elements. 

This evolution scenario was confirmed by the study of 
iHagenaail ()2001[ ). who analyzed a large statistical sam- 
ple of ERs using time series of magnetograms taken with 
the Michelson Doppler Imager (MDI) aboard the SOHO 
satellite. She found an average lifetime of ~ 4 hours, an 
average separation of ~ 9 Mm, depending on the total 
amount of magnetic flux carried by the ER, and a sepa- 
ration velocity decreasing from 4 km s~^ to 1.5 km s~^. 
The flux concentrations fragmented and moved apart, 
merging with other concentrations or with intranetwork 
flelds. 

The origin of ERs is quite controversial. They 
could represent the small-scale tail of the distribu- 
tion of activ e regions, gen erated by the l arge scale so- 
lar dynamo (|Parkedll955h . as argued by iHarvev et al.l 
( 19751). or be generated by a near- s urface local dy namo 
(jNordlund et"S][l99l : [PaAeJ[l99a [Cattaned[l99^ Al- 
ternatively, they might be the end product of failed active 
region e mergence processes that end in a catastrophic ex- 
plosion (More no Insertis et al.l[T995l). 

To shed light into their origin, iHagenaad (pOOl also 
analyzed the relationship of ERs with the solar mag- 
netic activity, concluding that only 60% of the ERs 
had an orientat i on co nsistent with Hale's polarity law. 
iHagenaar et al.l ()2003[ ) found that the number of ERs 
varies by a factor of 1.5 in antip hase with the solar cycle , 
confirming the early findings of lMartin fc Harveyl (|1979t ) , 
and that the latitude distribution of ERs is broader than 
that of the acti vity belts. The frequ ency spectrum of 
ERs deduced bv lHagenaar et all ()2003D appeared to be a 
continuous, smoothly decreasing distribution of bipolar 
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flux regions, formed by two distinct but coupled power- 
law spectra, with a turnover at around 2 x 10^^ Mx. This 
suggested the coexistence of both dynamo mechanisms, 
the global dynamo producing cycle-modulated large ERs 
and active regions, and the local dynamo generating 
small-scal e bipoles independent of the cycle phase. Fur- 
thermore, iHagenaar et al.l ()2008( ) showed that the rate of 
ERs emergence depends on the local flux imbalance: it 
is lower within strongly unipolar regions by a factor of 
3 relative to flux-balanced areas in the quiet Sun. This 
indicated that a large background field may affect the 
ERs generation process in the sub-surface layers. 

In the meantime, our knowledge of small-scale 
magnetic fields greatly improved from below the ERs 
threshold at 10^® Mx down to 10^^ Mx, at the limit 
of current detection capabilities. Thanks to the high- 
resolution observations taken by both ground-based 
telescopes with adaptive optics, and the Solar Optical 
Telescop e fSOT:lTsuneta et al.l l2008l) aboard the Hinode 
satellite (jKosugi et al.ll2Q07D . a number of small low- flux 
bipolar structures have been studied in recent years. 
Several contributions reported on the emergence of 
horizontally inclined magnetic fields with strengths of 
a few hundred G in the quiet S u n, mo s tly at granular 
scale (iMartfnez Gonzalez et al.l 120071: iCent eno et al.l 
2007t lOrozco Suarez et al.l l2008t iLites et al. 2008; 
Ishikawa et al.l 20081: Martinez Gonzalez fc Bellot Rubig 
20091: IGomorv et al.l I2010t Ishikaw aeTal.l 120101: 
Danilovic et al.l I2011D . The usual evolution of these 
structures begins with the appearance of an isolated 
patch of linear polarization signals, indicating horizontal 
fields, followed by the appearance of opposite-polarity 
circular polarization signals at opposite edges of the 
horizontal patch. These results are often interpreted 
as the rise of small-scale magnetic loops through the 
photosphere. 

Different mechanisms have been proposed to explain 
the decay of the bipoles: submergence of il-loops driven 
by convective downflows, flux dispersal resulting from 
horizontal convective motions with the Stokes signals 
falling below the instrumental sensitivity, ohmic dissi- 
pation, and the rise of U-loops towards the upper pho- 
tospheric layers following reconnection below the solar 
s urface. 

iParnell et cill (j2009[) studied MDI and SOT data to 
determine the flux distribution of all magnetic features 
that can be observed with these instruments on the solar 
surface at any one time. They found that flux concentra- 
tions can be described by a single power-law distribution 
with a slope of -1.85 for fluxes between 2 x 10^'' Mx to 10^^ 
Mx. These authors argue that the differences in flux dis- 
tributions found in previous works result from counting 
and identifying the magnetic structures in different ways, 
i.e., from selection effects. Their findings support a scale 
invariance of the distribution of fluxes that conceivably 
may extend to even smaller magnetic features beyond 
the current detection limit. Such a distribution might be 
produced either (i) by a solar dynamo that operates in 
the same way at a continuum of scales, and whose length 
scale increases with depth or (ii) if the magnetic features 
are generated by distinct dynamos with different charac- 
teristic lengths, by reprocessing of flux at the surface. In 
the latter case, such a process, through merging, cancel- 
lation, and fragmentation, would dominate to produce a 



p ower-law distribution . 

iThornton fc ParneUl (|201lD investigated the distribu- 
tion of small-scale intranetwork concentrations count- 
ing only the features related to emergence events with 
fluxes smaller than 10"'^^ Mx. By combining their re- 
sults with the results on the distributions of the emer- 
gence at larger scale, they found a single power law for 
all emerging features with fluxes ranging from 10^^ Mx 
to 10^'^ Mx wit h a slope of -2.7, which is consistent with 
the findings of Hagena ar et aP (|2003[ ). The emergence 
rate is dominated b y small-scale fields, in agreement with 
iSocas-Navarro fc Sanchez Almeidal (|2002l ). Since the sin- 
gle power-law distribution is already found at emergence, 
before the surface flux reprocessing, this result is sug- 
gestive of the idea that dynamo action occurs on all 
scales. Large magnetic structures would be created in 
the tachocline, whilst the smaller features would con- 
tinuously occur by a turbulent dynamo action over a 
range of scales throughout the convection zone. Nu- 
merical conyection simulations also support this scenario 
(|Stein et al.ll2010[) . The steeper slope of emerging struc- 
tures (-2.7) with respect to the magnetic features present 
at any instant on the solar surface (-1.85) can be ex- 
plained taking into account the shorter lifetime of the 
small-scale fields and their coalescence at the edges of 
the supergranules. 

In this paper, we report on the emergence of a magnetic 
bipole, observed at high spatial res olution by the Imaging 
Magn etograph experiment (IMaX; IMartfnez Fillet et al.l 
120 llh mounted on the 1-m aperture telescope flown 
on board the Sunrise balloon-borne solar observatory 
(iSolank i et al.H2010l: 'B arthol et al.ll2011l: iBerkefeld et al.l 
120111 : iGandorfer et al. 2 0111 ). We present a high-cadence 
analysis of the temporal evolution of this magnetic struc- 
ture. The flux content of the bipole is ^ 6 x 10^^ Mx, 
which places it close to the limit between the ERs and 
the small-scale fields observed at granular scale. In Sect. 
2 we describe the observations, in Sect. 3 we report our 
results on both the emergence and the decay phases, and 
in Sect. 4 we discuss our findings in a general context. 

2. OBSERVATIONS AND DATA ANALYSIS 

We have analyzed an IMaX data set obtained on 2009 
June 9, 01:30:54-02:02:29 UT. During that time, IMaX 
took polarization maps at five wavelength positions over 
the Fe I 525.02 nm line (Lande factor g = 3) a.t X = -8, 
-4, +4, +8, and -1-22.7 pm from the line center, with a 
spectral resolution of 8.5 pm. Six images were accumu- 
lated per wavelength point, and the full Stokes vector 
was recorded (V5-6 mode) . The temporal cadence of the 
scans is 33 s, with a pixel size of 0'.'055. Figure [T] shows 
the field of view (FoV) covered by the observations, about 
50" x 50" over a quiet region at the disk center. 

All data have been corrected for instrumental effects, 
by performing dark-current subtraction, flat-field correc- 
tion, and cross-talk removal. The blueshift over the FoV 
due to the coUimated setup of the Fabry-Perot etalon 
of the magnetograph is corrected in the inferred veloc- 
ity values. Two different types of data are produced: 
non-reconstructed data (level 1) and reconstructed data 
(level 2), obtained by using phase-diversity information. 
The reconstruction requires an apodization that effec- 
tively reduces the IMaX FoV down to about 45" x 45". 
In this study, we analyze level 1 data, since the emerging 
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Figure 1. Map of the non-reconstructed mean circular polariza- 
tion signal Va, i.e., the LOS magnetic flux, covering the full FoV 
of IMaX of about 50" X 50", obtained during the observations of 
the emerging bipole analyzed in this work. The square (dashed 
line), with a FoV of 10"5 X 10"5, indicates the sub-region whose 
spectra have been inverted with the SIR code. The inner rectangle 
(solid line) , with a FoV of ~ 4" X 8" indicates the location of the 
emerging flux region, as shown in the sequences displayed in Fig.|3] 



bipole was observed just at the right border of the full 
FoV (within the box in Fig. [T]). The spatial resolution is 
O'.'S (before reconstruction) and the noise level is about 
1 X 10~^ in units of continuum intensity per wavelength 
point in each Stokes parameter. Further details about 
data reduction are provided by iMartmez Fillet et al.l 

(Eon). 

We obtained maps of the mean circular polarization 
averaged over the line, Vg , and of the mean linear polar- 
ization signal I/s, given respectively by 
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and 



4 {Ic) ^ 



where (Ic) is the continuum intensity averaged over the 
IMaX FoV, e = [1,1,-1,-1], and i runs over the first 
four wavelength positions. 

We carried out inversions of the observed Fe I 
525.02 nm Stokes vector spectra u sing the SIR code 
(jRuiz Cobo &: del Toro Iniestalll992l ). This code numer- 
ically solves the radiative transfer equation along the line 
of sight (LOS) for Zeeman-polarized radiation under the 
assumption of local thermodynamic equilibrium. 

The inversion yields the temperature stratification in 
the range —4.0 < log r < 0, where r is the optical depth 
of the continuum at 500 nm. We use the Harvard Smith- 
sonian Reference Atmosphere (HSRA; iGingerich et al.l 
l(1971f) as initial model. The temperature is modified 
with two nodes. SIR also provides the LOS velocity 
vlos, the micro-turbulent velocity Vmicro-, the magnetic 



field strength B, the inclination and azimuth angles 7 
and (j) in the LOS reference frame. As these IMaX ob- 
servations were taken at disk center, the returned mag- 
netic parameters do not need to be converted to local 
solar coordinates. The magnetic filling factor has been 
assumed unity for these inversions, not taking into ac- 
count stray light contamination. The synthetic profiles 
are convolved with the spectral Point Spread Function 
at the focal plane of IMaX. 

In a first series of inversions, all the physical parame- 
ters are assumed to be constant with height, except for 
the temperature that has two nodes, as already men- 
tioned above. In another series of inversions, vlos and 
B are also allowed to have a linear dependence with the 
optical depth t. Such a gradient in both parameters has 
been found to be necessary in order to fit some profiles 
that show a considerable asymmetry between the blue 
and red lobes of Stokes V, or between the lobes of Q and 
U when they are significant. Elsewhere the fits obtained 
without gradient in velocity and magnetic field strength 
is found to be as good as the one with gradients. The 
total number of free parameters is 7 for SIR inversions (9 
if gradients in v^os and B are included) and the number 
of the data points in the V5-6 mode is 20. 

The absolute velocity is calibrated using the results 
coming from the Gaussian fits ob tained for the same 
data set (see, e.g.. lRoth et al.|[2010D . compared with the 
results of the SIR inversion without gradients. This cal- 
ibration uses a granular mean blueshift of 200 m s~^. 

In Fig. [2] we display a subset of the typical profiles 
acquired by IMaX, and the corresponding fits obtained 
using the SIR code. The spectra plotted in blue, show- 
ing predominantly horizontal fields, and red, with more 
vertical fields, show profiles well fitted without any gra- 
dient. There are also clear examples of spectra showing 
noticeable asymmetries in Stokes V in absence of a sig- 
nificant Ls signal (plotted green), as well as in Stokes 
Q and V with clear horizontal fields (plotted orange), 
which both require a gradient in v^os and B. These 
profiles show either downflows or upflows, respectively, 
as can be deduced also from the shifts with respect to 
the line center in Stokes /. The residuals of the fit, de- 
fined as Pobs — Pfit , with P being any of the polarization 
profiles, show an rms variation over the inverted FoV of 
1 — 1.6 X 10"'^, in agreement with the estimated noise per 
wavelength sample of the data. 

3. RESULTS 
3.1. Flux emergence phase 

Figure [3] shows a temporal sequence of the continuum 
intensity /c, and of the polarization maps and Vg. 
This sequence tracks the emergence of the bipole since 
its first detection in the photosphere. Green contours 
enclose the region with a significant linear polarization 
signal, i.e. Lg > 0.5%, whereas red and blue contours 
mark areas with positive and negative circular polariza- 
tion, respectively, with Vg > 0.8%. We also provide a 
movie of the evolution of the bipole in the online jour- 
nal. The temporal cadence of each frame is 33 s. 

At 01:35:54 UT the bipole appears in frame 2 at coor- 
dinates [3",3'.'5], inside a pre-existing granule visible at 
the same location in frame 1. It is later recognizable as 
a patch with increasing linear polarization signals in the 
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Figure 2. Observed (symbols) /, Q, U, and V profiles, for a sample of IMaX spectra, with the corresponding fits (lines) obtained with 
the SIR code. Blue (triangles/long dashed): horizontal fields, no gradients. Red (squares/dashed): vertical fields, no gradients. Green 
(circle/dot-dot-dot dashed): asymmetry in Stokes V with weak Ls signal, with gradients. Orange (diamonds/dot dashed): asymmetry 
in the blue and red lobes of Stokes U and V, with gradients. The dotted line represents the fit without gradients to the plotted orange 
spectrum (for Stokes / the dotted line falls on top of the orange one). The dashed vertical line indicates the nominal line center. 



following Ls maps that always coincides with the top of 
the granule during the whole emergence process. There 
is no evidence of any influence of the emerging magnetic 
flux on the granulation pattern. 

Note that the bipole emerges in a region where two 
subarcsecond bipoles were already present as seen in the 
first Vs frame of Fig. [3] (top panel). One of these bipoles 
is located at [3", 4"], nearly exactly over the region of 
emergence of the negative footpoint of our bipole, but 
with its polarities flipped (i.e., with the positive polarity 
toward the top part of the figure) with respect to those of 
the emerging bipole. The other is located at the bottom 
right of the studied area, at [3'.'5,2"], and possesses the 
same orientation as our bipole. 

The corresponding Vs footpoints of the large emerging 
bipole are first seen in the frames 3 - 4, as two small 
patches of opposite polarities at symmetrically opposite 
edges of the enhanced Ls patches. The evolution of the 
emerging loop can be easily followed in the subsequent Vs 
frames, where the opposite polarities of the bipole, i.e., 
the footpoints of the loop, are seen to separate from each 
other in opposite directions, with quite a strong Ls signal 
in between them being maintained until frame 9. This 
implies horizontal field in mid-photospheric layers, which 
is consistent with continuing flux emergence during this 
whole period of time. 

It is clear that the emerging bipole studied in this work 



interacts with the two pre-existing ones. The first bipole 
at [3", 4"] is completely overtaken by the negative polar- 
ity footpoint of our larger emerging bipole. The positive 
polarity footpoint of the new bipole, however, first en- 
ters an intergranular region but crosses over it to interact 
with the negative polarity footpoint of the pre-existing 
bipole (frame 5 in Fig. [5] at [3'.'5,2'.'5]). The negative 
pre-existing footpoint is completely washed away, while 
the emerging positive footpoint appears to aggregate the 
positive polarity flux of the pre-existing bipole. 

We plot the distance between the footpoints as a func- 
tion of time in Fig. S] {left panel), calculated as the 
distance between the centroids of the areas where \Vs 
is greater than 0.5%. The plot displayed in Fig. |4] 
{left panel) shows that the distance increases linearly 
during the first 1000 s of the loop evolution, with a 
mean separation velocity of 3.5 km s^^. Maximum foot- 
point separation is 4.5 Mm. This result has to be com- 
pa red with the analysis of emerging b i poles carried out 
bv lMartmez Gonzalez fc Bellot Rubiol (|2009| ). In some 
cases, they found an initial velocity separation of the 
footpoints of 6 km s~^, with a change after 500 s to w 2 
km s~^. Note that our bipole is comparable to the largest 
bipoles studied in that work. 

The corresponding temporal sequences of the physical 
parameters retrieved by SIR during the emergence are 
shown in Fig. [SJ for frames 1 - 8 of Fig. [3l and in Fig. IH 
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Figure 3. Top rows: Temporal sequence of the continuum maps during the emergence phase of the bipole. Time runs from left to right 
and continues in the lower set of panels. Red (blue) contours over the maps represent a circular polarization signal of +0.8 (-0.8) % of the 
Ic, and green contours represent a linear polarization signal of 0.5 % of the Ic- Middle rows: the same for the mean linear polarization 
signal Ls . Bottom rows: the same for the total circular polarization Vs ■ The elapsed time between consecutive frames is 33 s. From frame 
11 the decay phase begins. An animation of this figure is available in the online journal. 



6 



Guglielmino et al. 



for frames 9 - 16. We display the values of LOS veloc- 
ity and of the magnetic field strength averaged between 
log r = —1 and log r = —2, where the response functions 
are more sensitive (note that no filling factor is included 
in the inversions). Blueshifts, which indicate upward mo- 
tions, are present in the region since the beginning of the 
emergence until frame 10. They are cospatial with the 
region with strong Lg signals, indicative of nearby hori- 
zontal fields, and correspond to the granule in the con- 
tinuum map. The emerging flux region is surrounded by 
downflows along the intergranular lanes which outline the 
granule. The larger scale variable velocity pattern near 
the top of the velocity frames 5 - 7 reflects t he p- mode 
oscillations, as discussed by, e.g.. lRoth et al.l ()201ClD . 

The magnetic field strength has a peak of about 400 
G around frames 5-6, then it begins to fade (note that 
this occurs before the footpoints enter the intergranular 
lanes, which happens in frame 7). As our filling factor 
is unity, this value should be taken as a flux-equivalent 
field strength. The maximum values of B are found in the 
region where the magnetic field is more horizontal. These 
strengths are of the order of the typical equipartition field 
strength Bg for granules, given by 

Be^ 1 

where /io is the magnetic permeability of free space, with 
t> pa 2 km s^^ and p — "i x. kg m~^ (jlshikawa et al.l 
[20081) . 

The footpoints do not reach a vertical orientation, but 
they remain rather inclined at an angle of about 45° with 
respect to the vertical. Furthermore, they surround the 
granule having a crescent-shaped aspect, different from 
the usual circular shape found in previous observations. 

In the emergence zone, the azimuth angle is quite ho- 
mogeneous at every instant during the rise, with the scat- 
ter ranging from ±5° to ±15° around the average value. 
However, the mean value, indicated by arrows in Figs. [5] 
and ini changes with time, rotating counterclockwise 10° 
per frame until the end of the emergenc^l. The observed 
rotation rate is thus 0°.3 s""'^. In the Lg and Vg maps in 
Fig. [21 we can see that the bipole axis also rotates in the 
same counterclockwise direction. This similar rotation 
for both the field azimuth angle and the bipole footpoints 
is compatible with the absence of writhe in the emerging 
flux tube. The observed evolution of the azimuth, which 
steadily varies by IS.(j) « 90° during the emergence, may 
be attributed to convective motions, but we cannot rule 
out the presence of some twist in the structure, unfolding 
during its ascent. 

The total magnetic fiux content is given by 

$ = ^B,cos7jS' , 

i 

where index i runs over all the pixels of the bipole with 
the same polarity and 5* is the area of the Sun covered by 
a pixel. We note here that our assumption of a filling fac- 
tor equal unity, while probably affects the values found 

* Note here that we do not solve the azimuth ambiguity. The 
arrow is selected simply by choosing the azimuth direction that 
points from the positive to the negative footpoint of the bipole 



for the magnetic field strength, has no impact on the es- 
timate of the magnetic flux density and on the total flux 
content. In both polarities, the flux grows linearly until 
a maximum of about 6 x 10^^ Mx at frame 8 is reached, 
4 minutes after the first detection of the bipole as seen 
in Fig. [4| {right panel). In this figure, the fluxes were 
estimated by adding all the contributions from manually 
selected pixels in the footpoint area of interest. Miss- 
ing points correspond to non-inverted frames. The flux 
growth rate is thus 2.6 x 10^^ Mx s~^. 

Note that the slightly smaller positive fiux (red sym- 
bols) may be due to a different inclination between the 
footpoints and, as footpoints move apart, to the fact that 
the corresponding footpoint is located closer to the right 
boundary of the observed area. At the end of the emer- 
gence process (defined as the maximum in the flux curve) , 
the distance between the footpoints of the bipole is about 
1.5 Mm, with roughly the same size as a typical granule. 
All these results are consistent with the emergence of an 
ri-shaped magnetic bipole at granular scale. 

The maximum distance D attained by the foot- 
points, 4.5 Mm, is roughly in agreement with the linear 
log -D/ log $ relation found by Hagenaar (2001) for ERs, 
which gives D ~ 3.5 x ^^-^^ = 5.5 Mm. 

Active region flux emergence displays upfiows as- 
sociated with transverse fields at equipartition values 
and downfiows at the footp oints with much larger field 
strengths (see, e.g.. lLites et al. 1998; Solanki ct al. 2003). 
In order to compare this behaviour with our bipolar 
emergence, we have produced a scatter plot of the mag- 
netic field strength and LOS velocity vs. zenith angle 7, 
as shown in Fig. [7| for the SIR inversion without gradi- 
ents in B and v^os (frame 12). In that frame, the foot- 
points of the loop, which correspond to the regions with 
more vertical field, have also stronger field strengths of 
about 200 — 250 G with downflows slightly larger than 1 
km s""'^, decreasing with time. As already stated above, 
it is conceivable that these fields are indeed at equiparti- 
tion values. In the horizontal field regions, an upflow of 
1.0 km s^^ is seen. This upward motion also decreases 
with time. Peak upflows for the transverse fields are seen 
in frames 7-9 where they reach 2.5 km . 

No evidence of asymmetry between the footpoints of 
the bipole is found in the data in both the field strength 
and the LOS velocity. However, the linear polarization 
signal, Lg is concentrated mainly at the negative foot- 
point in frames 10 - 13. 

3.2. Flux decay phase 

The observations used in this work allow a detailed 
analysis of the decay phase of the bipole, the less known 
stage of the evolution of small bipoles and ephemeral 
regions. The footpoints of the emerged bipole remain 
visible until the end of the emergence process and well 
beyond into the decay phase. This phase can be defined 
to start at the time when the total measured flux be- 
gins to decrease, i.e., from frame 9 onwards. Note that 
the flux curves in Fig. [4| {right panel) never reach a flat 
portion where the flux of the footpoints would stay con- 
stant. Instead, as soon as the maximum is reached, a 
mechanism somehow begins to erode the newly emerged 
flux. 

Figure [H which is the continuation of Fig. [3j shows 
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Figure 4. Left panel: Distance between tlie centroids of the opposite polarity footpoint of tlie emerging bipole as a function of time. Tlie 
dotted vertical line indicates the time at which separation deviates from a linear increase. Right panel: Evolution of the flux content in 
the positive (negative) polarity of the emerging bipole. Red (blue) symbols indicate data values, solid lines represent offset exponential fits 
to data during the decay phase. 



77 5 

CD 

ffi A 
U ^ 

^ 3 
2 

1 





1 1 1 1 1 1 1 

LOS; veloc 

" 1. , . , . ., . ' 


1 1 1 1 1 1 1 

2 , ^ , , , / 


, 1 1 1 .1 1 , 

c 



..^.■/i..u., 


1 1 1 1 1 1 1 

4 


.1 1 

■ • ' 


1 ' 1 ' 

6',', , 


1 1 1 1 .1 1 1 


1 1 , 1 1 1 1 

















7 Inclinotion 




D 



t3 















2 J 
[arcsecl 



: 2 3 

[arcsecl 



1 2 .3 
[arcsecl 



2 .3 
[arcsecl 



z o u 



[ 



-0 



I 



r4oo 

-300 

cn 

CO 

-200 I 
o 

- 100 

.0 
180° 
135° 
90° 

45° 
0° 

180° 
-90° 
-0° 



i 



-90° 
-180° 



Figure 5. LOS velocity and magnetic field strength at log r = —2 (first and second rows), inclination and azimuth angles (third and 
fourth rows), for frames 1-8 (compare with Fig.lSj. Contours represent a linear polarization signal of 0.5 % of the Ic, as in Fig.|3] Arrows 
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the evolution of the bipole after the emergence process: since the evolution during this phase is less rapid, only 
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Figure 6. Same as in Fig. [5] for frames 9 - 16. 

every other frame is shown (but see also the on-hne ma- 
terial). Footpoints move apart until they start breaking 
into smaller flux patches and interact with the surround- 
ing fields. Their remnants remain distinguishable at least 
for 15 minutes after the first detection of the magnetic 
structure. In frames 16 - 18, the bulk of the negative 
polarity starts being fragmented, losing its compactness. 
During this process the various patches of this polarity 
move along intergranular lanes. 

In frames 32 - 34 the negative polarity of the bipole 
totally merges into a pre-existing flux element of the same 
polarity, and is no longer distinguishable as an individual 
flux concentration. This explains why the flux curve for 
the negative polarity footpoint in the right panel of Fig. 2] 
stops there. 

The positive polarity footpoint at the beginning of the 
decay phase has merged with the pre-existing positive 
flux of the bipole at the bottom. Note that, in this case, 
no net flux is lost as the positive polarity footpoint of the 
pre-existing bipole is aggregated to the positive polarity 
of the newly emerged one. This could explain why this 
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cancellation episode is not apparent in the red flux curve 
of Fig. m (ri^/if panel). After this cancellation, the newly 
created positive polarity patch is seen to move near the 
edge of our field of view with little evident interaction 
with other fluxes (see Fig. [8]) . Interestingly, during this 
time the positive fiux patch moves largely over granula- 
tion dominated regions instead of intergranular ones. 

In order to get an indication of the physical processes 
occurring during the decay phase, we have fitted the fiux 
evolution curves during the decay period using the MP- 
FIT routine (Markwardt 2009), which performs fits to 
data with user-defined functions. We have used as a 
guess function an exponential of the type 

f{t)=a + b-e'^* . 

We obtained ;3± = -2.5 xlO"^ and-S.OxlQ-^ s~i 
for the positive and negative fiux, respectively. This cor- 
responds to a "turbulent" diffusion costant rj — L'^/t of 
the order of 8 x 10^ km^ s~^, where r = /3~^ is the e- 
folding time and L « 0'.'8 is the size of the footpoint at 
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Figure 7. Scatter plots of B and fj^oSi ^ ^ function of tiie inclination 7, for a representative frame (No. 12 in Fig. [Sjl. We show the 
results for the SIR inversion without gradients, and the corresponding maps of the physical parameters. Only the pixels within the white 
contour in the B map with a signal than three times the noise have been indicated in the plot. Keep in mind that inversion have a tendency 
to retrieve horizontal fields for pixels at noise level. 



the time of the flux peak in Fig. 2] {right panel). It is 
interesting to note that this diffusion constant compares 
rather well with those inferred by flux transport models 
of the large scale field, typica lly set at 6 x 10^ km^ s~^ 
(see the instructive review bv lSheelevI [20051) . 

This order of magnitude agreement between the diffu- 
sion constants of the large scale field and of our granular 
scale bipole is rather illuminating. As the flux loss oc- 
curs in situ, at a rate similar to that encountered in flux 
transport models, this could indicate an analogous phys- 
ical process acting in both flux loss scenarios. 

Even though the flux decay curve is clearly non-linear, 
it is interesting to compute a mean flux decay rate which 
turns out to be 2.86 x lO^"* Mx s^^. This rate is one 
order of magnitude smaller than the flux growth rate. A 
similar asymmetry between the gr owth and decay rates is 
known to e xist for active regions (|Martmez Fillet 1 120021 : 
Ivan Driel-G csztelvi 2002). 

4. DISCUSSION AND CONCLUSIONS 

IMaX/Sunrise observations of the solar photosphere 
taken at disk center have revealed a number of small- 
scale episodes of niagnetic flux emerge nce (see, e.g., 
iDanilovic et al.ir2010t iSolanki et al.|[20Tol and references 
therein). We have carried out an analysis of the emer- 



gence and disappearance of a small magnetic bipole. We 
have analyzed the polarization maps and then we have 
inverted the Stokes profiles with the SIR code, to obtain 
information on the physical parameters of this magnetic 
structure. 

The detection of weak field strengths (400 G) at 
emergence with horizontal inclination, associated with 
blueshifted Stokes line profiles between opposite polarity 
Stokes V profiles, is a well-known signature of the emer- 
ge nce of a concentra ted flux loop, as firstly pointed out 
bv lLites et al.l (jl996D at active region scales. At smaller 
(granular) scales, our analysis shows results resembling 
those found by Gomory et al. (2010), but with a higher 
cadence and better spatial resolution. 

The magnetic flux content places this small bipole at a 
halfway p o int be tween the ephemeral regions studied by 
iHagenaan ()2001| ) and the all-pervasive loops of the quiet 
Sun. The bipole indeed appears to be coincident with 
a granule in the continuum map, so it would apparently 
represent a typical case of flux emergence at granular 
scale. However, the present case is the largest structure 
observed to emerge in the IMaX/Sunrise data analyzed 
so far. 

Once emerged, the magnetic fiux spreads over an area 
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Figure 8. Same as in Fig. \3\ but during the decay phase of the bipole. Frame numbering continues from the last frame of Fig. [S] In this 
figure, the elapsed time between consecutive frames is 66 s. 
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larger than the source granule. The footpoints separate 
up to 4.5 Mm from each other and their remnants remain 
visible until 15-20 minutes after the first detection of 
the magnetic structure. Nevertheless, the granulation 
pattern seems not to be affected by the flux brought to 
the surface by the emerging b i pole. This agrees with the 
simulations of iCheung et al.l |2007[) . showing that flux 
tubes carrying less than ~ 10"'^^ Mx of longitudinal flux 
do not produce visible disturbances to the granulation 
pattern. 

During the emergence, the footpoints of the bipole do 
not have the more t ypical circular shape seen in Hinode 
observations (e.g., iMartmez Gonzalez fc Bellot Rubiol 
|2009( ). They rather surround the granule and have a 
crescent-shaped appearance. Moreover, footpoints do 
not reach a vertical orientation: they form an angle of 
about 45°. This could suggest that only a part of the 
O-loop is observ ed above the continuum formation layer, 
as also found bv llshikawa et al.l ()2010( ). or that the field 
lines follow the borders of the granule, which could not 
be vertically steep. Another possibility would be that the 
field in the footpoint of this emerging loop is not resolved 
at the spatial resolution ot IMaX level 1 data. Former 
measurements presented bv lSolanki et al.l ()1996f ) showed 
that quiet-Sun feature with an amount of flux over 10^'' 
Mx tend to have kG field strength, while the inversions 
retrieve a field strength of « 250 G. An underestimate 
of the field stren gth im plies that the inclination angle 7 
is overestimated (jStenfl o 1985a, bi), so that the magnetic 
field in the footpoints would be more vertical than the 
inversions indicate. Anyway, such a difference of more 
than a factor of 3 in the field strength is rather diffi- 
cult to conceive for our measurements based on the high 
Zeeman sensitive Fe I 525.02 nm line (g = 3) displaying 
Stokes Q and U signals above the noise level (see, e.g., 
frame 12 in Fig. |3]). These signals help inferring the field 
strength in the inversions more reliably than when only 
Stokes V is used. 

The azimuth angle changes during the emergence by 
about 90° (see Figs. [S] and [5]), in the same direction as the 
line joining the footpoints. Convective motions acting on 
the emerging loop or the unfolding of any twist present 
in the original structure could generate this behaviour. 

The fiux decay follows an exponential law. The 
observed cancellation and merging with a pre-existing 
bipole does not have an effect on the flux history of the 
bipole and is not the dominant decay process in this ex- 
ample, although it might change the connectivity of the 
flux patches. Flux removal dominates from the moment 
in which the maximum flux is observed. Clearly, flux is 
destroyed in situ, i.e., without visible interactions with 
opposite polarity patches, in some efficient way that dom- 
inates the fiux budget as soon as no more fiux is brought 
to the surface. This elusive process probably acts also 
during the fiux emergence phase, but is not evident at 
this stage simply because emergence is one order of mag- 
nitude stronger than the flux decay rate. The diffusion 
constant estimated for this process is only a little bit 
higher than typical estimates used for large scale flux 
transp ort models (6 x 10^ km^ I Cameron et al.l 

(|2011l ) have also found similar values in numerical sim- 
ulations of small-scale mixed polarity field in the near- 
surface layers of the Sun. They show that magnetic el- 



ements are advected by the horizontal granular motions 
against each other and that the flux is removed at a 
rate corresponding to an effective turbulent diffusivity 
of about 3 X 10^ km^ s~^, which is only slightly smaller 
than our estimate. 

The existence of flux-removing diffusive processes that 
are similarly effective at active region and granular scales 
deserves to be further studied from a theoretical and ob- 
servational point of view. Our value of the diffusion con- 
stant can be used to esti mate the time in wh ich an active 
region would decay. As iMever et al.l ()1974[ ) inferred for 
the turbulent diffusion of a monolithic magnetic flux con- 
centration, the decay time is given by 



t 



where is the original flux content, i] the turbulent dif- 
fusivity, and Be the equipartition field strength. Using 
this formula for a typical active region of 10^^ Mx with 
our value of the turbulent diffusivity (8 x 10^ km^ 
we estimate a decay time of « 3 days. As active regions 
live longer, it seems clear that the diffusivity found in 
this paper is more effective on smaller flux structures, 
or possibly on flux structures embedded in the vigorous 
granular convection found in the quiet Sun, rather than 
in the smaller and more slowly evolving abnormal gran- 
ulation in active regions. 

Finally, it is interesting to compare the emergence rate 
with the distribution of e merging fluxes proposed by 
[Thornton fc Parnelll (|2011f) . The number of events per 
unit area per unit time within a given flux range $ is 



iVev ($) 



no 



a V <i>o 



$1 



where no = 3.14 x 10"" cm'^ day"^ and $0 = 10^^ Mx. 
Remembering that they found a power index a = —2.7, 
we obtain a frequency of A^ev = 4.9 x 10^^^ cm^'^ day~^ 
for a flux range from 5.5 x 10^'' to 1 x 10'^^ Mx. Taking 
into account our FoV of ~ 1300 Mm^ and the duration of 
the observations of 30 minutes, we get A'ov — 14 events 
(i.e., 7 bipolcs). The observed number of flux patches 
within this rang e in our data set is smaller t han the pre- 
diction from the lThornton fc Parnelll (j201lD work, as we 
found only the two footpoints of the bipole. Note also 
that no flux patches with fluxes above 6.5 x 10^^ Mx are 
observed, whereas the above distribution would have sug- 
gested otherwise. Given the short duration of the period 
analysed here, it is unclear whether this disagreement is 
fortuitous or indicates an excess in the predicted number 
of bipoles due to an overestimate of the total density of 
emergence events. The activity minimum in which the 
Sun was residing during the Sunrise flight could also 
be behind our low number of observed bipoles. Note 
that the difference in both time and space resolution be- 
tween the two data sets could affect this comparison: 
th us, it would be worthwhile a pplying the algorithm used 
bv lThornton fc Parnelll (|2011D to these Sunrise observa- 
tions. 

On the other hand, the bipole closely behaves as a 
small-scale version of an ER: we found a good agreement 
with the dependence of the footpoint separation with the 
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flux content proposed by iHagenaaii (|2001[ ). The overall 
morphology is quite similar, even if we do not find any ev- 
idence of mixed polarities in between the two footpoints 
during the emergence, as reported bv iGugli elmino et alJ 
(Horn) for an ER with flux content - 1.5 x lO^^ Mx, i.e., 
twenty times larger than the bipole studied here. 

Further studies of similar time series can be of great 
value to consolidate the results presented here. More 
cases are needed to clarify these findings, using spec- 
tropolarimetric data of a similar resolution, sensitiv- 
ity and temporal coverage such as those that could 
be obtained with specially tailored observations with 
SOT /Hinode, CRISP/SST, and IBIS/DST. We are not 
aware of such results yet and, thus, it is difficult for the 
moment to extend and compare these findings with more 
examples, and to understand their implications for the 
global solar flux budget. 
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